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ABSTRACT 



Prompted by a high optical state in September 2007, the Whole Earth Blazar 
Telescope (WEBT) consortium organized an intensive optical. ncar-IR (JHK) and 
radio observing campaign on the intermediate BL Lac object 3C 66A throughout 
the fall and winter of 2007 - 2008. In this paper, we present data from 28 
observatories in 12 countries, covering the observing season from late July 2007 
through February 2008. 

The source remained in a high optical state throughout the observing period 
and exhibited several bright flares on time scales of ~ 10 days. This included an 
exceptional outburst around September 15 - 20, 2007, reaching a peak brightness 
at ~ 13.4. Our campaign revealed microvariability with flux changes up to 
|di?/dt| ~ 0.02 mag/hr. 

Our observations do not reveal evidence for systematic spectral variability 
in the overall high state covered by our campaign, in agreement with previous 
results. In particular, we do not find evidence for spectral hysteresis in 3C 66A 
for which hints were found in an earlier campaign in a somewhat lower flux state. 
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We did also not find any evidence for spectral lags in the discrete correlation 
functions between different optical bands. 

We infer a value of the magnetic field in the emission region of i? ~ 
19 e^^'^ T^^^'^ dI^^"^ G, where is the magnetic field equipartition fraction, Th 
is the shortest observed variability time scale in units of hours, and Di is the 
Doppler factor in units of 10. From the lack of systematic spectral variability, we 
can derive an upper limit on the Doppler factor, D < 28 t^^^^ ^^b^^ ■ This is in per- 
fect agreement with superluminal motion measurements with the VLBI/VLBA 
of Papp < 27 and argues against models with very high Lorentz factors of F > 50, 
required for a one-zone synchrotron-self-Compton interpretation of some high- 
frequency-peaked BL Lac objects detected at TeV 7-ray energies. 

Subject headings: galaxies: active — BL Lacertae objects: individual (3C 66A) 
— gamma-rays: theory — radiation mechanisms: non-thermal 



Introduction 



Blazars are the most violent class of active galactic nuclei, consisting of fiat-spectrum 
radio quasars (FSRQs) and BL Lac objects. They exhibit rapid variabili ty down to time 
scales as short as a few minutes (lAharonian et al.l 120071 : lAlbert et al.ll2007bl ). Their observed 
fiux is dominated by a non-thermal continuum exhibiting two broad spectral bumps: A 
low-frequency bump from radio to UV - X-ray frequencies, and a high-frequency component 
from X-ray to 7-rays. 

In the framework of relativistic jet models, the low-frequency (radio - optical/UV) 
emission from blazars is interpreted as synchrotron emission from nonthermal electrons in a 
relativistic jet. The high-frequency (X-ray - 7-ray) emission could either be produced via 
Compton upscattering of low frequency radiation by the same electrons responsible for the 



synchrotron emission (leptonic jet models; for a recent review see, e.g.. iBottcherl 120071 ). or 



due to hadronic processes initiated by relativistic p rotons co-accelerated with the electrons 
(hadr onic models, for a recent discussion see, e.g., 
2003h . 



Miicke fc Protheroe 2001: Miicke et al. 



The blazar 3C 66A (= 0219+428) is classified as a low-frequency peaked BL Lac object 
(LBL), a class also commonly referred to as radio selected BL Lac objects. Its nonthermal 



^The radio-to-optical data presented in this paper are stored in the WEBT archive; for questions regarding 
their availabihty, please contact the WEBT President Massimo Villata (villata@oato.inaf.it). 
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low-frequency spectral component extends from radio frequencies through soft X-rays and 
typically peaks in the optical frequency range. The high-frequency compo nent seems to 
peak i n the multi-MeV - GeV energy range. Since its optical identification by lWills &: Wills 
( I1974J ). 3C 66A has been the target of many radio, IR, optical, X-ray, and 7-ray observations 
in the past, although it is not as regularly monitored at radio frequencies as many other 
blazars due to problems with source confusion with the nearby rad io galaxy 3C 66B (6'. 5 
from 3C 66A), in partic ular at lower (4.8 and 8 GHz) frequencies (lAUer. AUer. fc Hughes 



1994; Takalo et al. 



19961). 



To date, about 2 dozen blazars have been detected at very high energies (> 100 GeV) 
with ground-based air Cerenkov telescope facilities. Until very recently, all TeV blazars 
belonged to the sub-class of high-frequency peaked BL Lac objec ts (HBLs). However, the 



recent detections of the intermediate BL Lac obj ect W Comae (lAcciari et al.ll2008l ). the 
low-frequency peaked BL Lac object BL Lacertae (lAlbert et al.ll2007al ). and even the flat- 
spectrum radio quasar 3C 279 (lAlbert et al.ll2008l ) demonstrate the potential to extend the 
extragalactic TeV source list to all classes of blazars, in particular intermediate and low- 
frequency peaked BL Lac objects (LBLs) with lower z/Fj, peak frequencies in their broadband 
spectral energy distributions (SEDs). 

3C 66A has been suggested for quite some time as a promising candidate for detection 
by t he new generation of atmospher ic Cerenkov telescope facilities like MAGIC or VERITAS 
(e.g. ICostamante fc Ghisellinill2002l ) (it is too far north in the sky to be observed by HESS), 



and intensive observations of 3C 66A by the VERITAS collaboration have recently resulted 
in the significant detectio n of VHE 7-ray emission from 3C 66A in September - October 
2008 JSwordv et al.ll2008h . 



In the optical, 3C 66A is generally observed as a point source, wi th no indication of the 
host g alaxy. The host galaxy of 3C 66A was marginally resolved by IWurtz. Stocke. fc Yee 
(119961 ). They found -Rcunn = 19.0™'^^ for the host galaxy; the Hubb le type could not be deter- 
mined . In 3C 66A, a weak Mg II emission line has been detected by lMiller. French. &: Hawley 
(119781 ). This led to the determination of its redshift at z = 0.444, which was later confirmed 



by th e detection of a weak Lya line in the lUE spectrum of 3C 66A (ILanzetta. Turnshek. &: Sandoval 
I993I ). However, as recently pointed out bv iBramel et al. (|2fl05[), these redshift determina- 



tions are actually still quite uncertain (see also iFinke et al.l l2008al ). In this paper, we do 
base our analysis on a redshift value of z = 0.444, but remind the reader that some results of 
the physical interpretation should be considered as tentative pending a more solid redshift 
determination. 



The long-term variability of 3C 66A at near-infrared (J, H, and K bands) and optica l (U, 
B, V, R, I) wavelengths has recently been compiled and analyzed by iFan Sz LinI (119991 ) and 
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Fan fc LinI (120001 ). respectively. Variability at those wavelengths is typically characterized by 
variations over <,1.5 mag on time scales ranging from ~ 1 week to several years. A positive 
cor relation between the B - R color (spectral hardness) and the R magnitude has been found 
by IVagnetti et al.l (|2003[). T he most recent multiwavelength campaign on 3C 66A in 2003 
- 2004 jBottcher et allboosh found hmts for spectral hysteresis, with the B - R hardness 
peaking several days prior to the R- and B-band fluxes during large fl ares. That ca r npaig n 
also confirmed the results from an intensive long- monitoring effort by iTakalo et al.l (119961 ) , 
revealing evidence f or rapid microvariability, including a decline ~ 0.2 mag within ~ 6 hr. 
Lainela et al.l (119991 ) also report a 65-day periodicity of the source in its optically bright 
state, which has so fa r not been confirmed in any other analysis. However, the campaign of 
Bottcher et al.l (120051 ) showed several major outbursts separated by ~ 50 - 57 days, possibly 
also indicating a quasi-periodic behavior. 

Superlurninal m otion of individual radio components of the jet has been detected by 
Jorstad et al.l ( 200111. While t he ide ntification of radio knots across different observing epochs 
is not unique, iJorstad et al.l (120011 ) favor an interpretation implying superluminal motions 
of up to /3app ~ 19 ^ 27, where h = iJo/(100 km s~^ Mpc~^) parameterizes the Hubble 
constant. This would imply a lower limit on the bulk Lorentz factor of the radio emitting 
regions of F > 27. Theoretical modeling using a time-dependent leptonic jet model produced 
an acceptable fit to the SEP and light curve o f 3C 66A during the 2003 - 2004 campaign 
with such a choice of F (jjoshi fc Bottcherll2007l ). 



In the fall of 2007, 3C 66A was found in a very active state, reaching a peak brightness 
around Sept. 14 of i? ~ 13.4. This triggered a new optical - IR - radio observing campaign 
by the Whole Earth Blazar Telescope (WEBli) collaboration with intensive monitoring 
throughout fall and winter 2007 - 2008. In this paper we present collected data from late July 
2007 through early March 2008. This high optical state also triggered very-high-energy 7-ray 
observations with the VERITAS array. Continued observations by VERITAS t hroughout the 



obser ving season in 2008 resulted in the VHE 7-ray detection mentioned above (ISwordy et al. 



20081 ). Details of those observations will be published in a separate paper. In the following, 
we will present in S|2] a summary of the observations and data analysis, and present light 
curves in the various frequency bands. In §3] we test for spectral variability to be derived 
from our results. Cross-correlations between the variability in different optical bands are 
presented in ^ We discuss our results and derive limits on the magnetic field and the 
Doppler factor in ^ and we summarize in ^ 



Throughout this paper, we refer to a as the energy spectral index, F^, [Jy] oc z/ A 



^ |http : //www ■ oato . inaf . it/blazaxs/webt/ 
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cosmology with = 0.3, = 0.7, and Hq = 70 km s ^ Mpc ^ is used. In this cosmology, 
and using the redshift of z = 0.444, the luminosity distance of 3C 66A is = 2.46 Gpc. 



2. Observations, data reduction, and light curves 

3C 66A was observed in a coordinated multiwavelength campaign at radio, near-IR, 
and optical frequencies by the WEBT collaboration during the 2007 - 2008 observing sea- 
son, from late July 2007 through early March 2008. The object is being continuously moni- 
tored at radio to o ptical wavelengths by the GLAST-AGILE Support Program (GASP; see 



Villata et al.ll2008l ). whose observers contributed also to this campaign. The overall time- 
line of the campaign, along with the measured long-term light curves at radio, infrared, 
and optical frequencies is illustrated in Fig. [H Table [1] lists all participating observatories 
which contributed data to this campaign. In this section, we will describe the individual 
observations in the various frequency ranges and outline the data reduction and analysis. 



2.1. Optical and infrared observations 



In the optical component of the WEBT campaign, 22 observatories in 12 countries con- 
tributed 3461 individual photometric data points. The observing strategy and data analysis 
followed to a large extent the standar d procedure used for the previous, successful WEBT 
campaign on 3C 66A in 2003 - 2004 ( iBottcher et al.ll2005l ). Observers were asked to per- 
form bias (and, where necessary, dark) corrections as well as flat-fielding on their frames, 
and obtain instrumental magnitudes, applying either aperture photometry (using IRAF or 
CCDPHOT) or Gaussian fitti ng for the source 3C 66A an d the comparison stars no. 13, 
14, 21, and 23 in the tables of I Gonzalez- Perez et al.l (120011 ) . where high-precision standard 
magnitudes for these stars have been published. This calibration has then been used to 
convert instrumental to standard photometric magnitudes for each data set. In the next 
step, unreliable data points (with large error bars at times when higher-quahty data points 
were available) were discarded. Our data did not provide evidence for significant variability 
on sub-hour time scales. Consequently, error bars on individual data sets could be further 
reduced by re-binning on time scales of 20 min. Finally, there may be systematic offsets 
between different instruments and telescopes. Wherever our data sets contained sufficient 
independent measurements to clearly identify such offsets, individual data sets were cor- 
rected by applying appropriate correction factors. The resulting offsets, by which individual 
data sets have been corrected to create a uniform light curve, are listed in Tab. [21 
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Fig. 1. — Timeline of the WEBT campaign on 3C 66A in 2007 - 2008, including the observed 
(before extinction correction) optical (UBVRI), near- infrared K-band, and radio light curves 
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Table 1. List of observatories that contributed data to this campaign. In addition to the 
standard optical Johnson's-Cousin's UBVRIJHK bands, wc obtained optical/near-IR data 
in the c, i, and o bands. The central wavelengths of these bands are c = 4206 A, i = 6658 

A, and o = 9173 A. 



Observatory Specifications frequency / filters / energy range N^ha 



Radio Observatories 



Crimean Radio Obs., Ukraine 


22 m 


36 GHz 


5 


Mcdicina, Italy 


32 m 


5, 8, 22 GHz 


19 


Metsahovi, Finland 


14 m 


37 GHz 


32 


Noto, Italy 


32 m 


43 GHz 


19 


UMRAO, Michigan, USA 


26 m 


4.8, 8, 14.5 GHz 


12 


Infrared Observatories 


Campo Imperatore, Italy 


1.1 m 


J, H, K 


72 


Optical Observatories 


ARIES, Nainital, India 


104 cm 


B, V, R, I 


15 


Armenzano, Italy 


35, 40 cm 


B, V, R, I 


70 


Belogradchik, Bulgaria 


60 cm 


V, R, I 


17 


Obs. de Bordeaux, Prance 


20 cm 


V 


34 


Catania, Italy 


91 cm 


U, B, V 


60 


Crimean Astr. Obs., Ukraine 


70 cm 


B, V, R, I 


147 


Gualba Obs., Spain 


35 cm 


V, R, I 


12 


Jakokoski Obs., Finland 


51 cm 


R 


69 


Kitt Peak (MDM), Arizona, USA 


130 cm 


U, B, V, R, I 


333 


L'AmpoUa, Spain 


36 cm 


R 


396 


Michael Adrian Obs., Germany 


120 cm 


R 


180 


New Mexico Skies Obs., USA 


30 cm 


V, R, I 


18 


Perugia, Italy 


40 cm 


R 


25 


Roque (KVA), Canary Islands 


35 cm 


R 


318 


Rozhcn, Bulgaria 


200 cm 


U, B, V, R, I 


139 


SabadcU, Spain 


50 cm 


V, R 


464 


San Pedro Martir, Mexico 


84 cm 


R 


10 


St. Petersburg, Russia 


38 cm 


B, V,. R, I 


50 


Torino, Italy 


105 cm 


B, V, R, I 


77 


Tuorla, Finland 


103 cm 


R 


19 


Valle d'Aosta, Italy 


81 cm 


B, V, R, I 


53 


Xinglong, China 


60 cm 


c, i, o 


699 
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Fig. 2. — Observed light curves in the optical R, c, i, and o bands (before extinction correc- 
tion) over the entire duration of the campaign. The central wavelengths of these bands are 
c = 4206 A, i = 6658 A, and o = 9173 A. 
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In order to provide information on the intrinsic broadband spectral shape (and, in partic- 
ular, a reliable extraction of B - R color indices), the data were de -reddened. For this p urpose, 
Galactic Extinction coefficients were calculated using Table 3 of ICardelli et al. Jl989h . based 



on Ab = 0.363 mag[§ and Ry = 3.1. As mentioned in the introduction, the R magnitude 
of the host galaxy of 3C 66A is ~ 19 mag, so its contribution is negligible compared to the 
average AGN magnitude of -R<,14, and no host-galaxy correction was applied. 

As a consequence of the chosen observing strategy, the R- and B-band light curves are 
the most densely sampled ones. The R-band light curve over the entire duration of the 
campaign is compared to the light curves at all other optical bands as well as the near-IR 
K-band and 36/37 GHz radio light curves in Figs. [1] and [21 These figures illustrate that the 
object remained in a very bright and active state, with an average magnitude of ~ 14.0, 
throughout the campaign. The source exhibited several large fiares with Ai?>^0.5 mag out 
to -R ~ 13.5 - 13.6, on characteristic time scales of ~ 10 days. The most dramatic outburst 
occurred aro und September 14, w ith a peak magnitude of i? ~ 13.4. As already found in 
2003 - 2004 jBottcher et al.lboosh . visual inspection of individual major outbursts suggests 



periods of more rapid decline than rise (see, e.g.. Fig. [3]). 

Our campaign revealed clear evidence for rapid intraday microvariability. Two examples 
are shown in Figures [3] and HI revealing a rapid decline of the optical brightness by dR/dt ~ 
0.016 mag/hr on JD 2454331 (August 18, 2007), and a rapid rise by dR/dt ~ -0.008 mag/hr 
on JD 2454402 (October 28, 2007). Our results do not provide any evidence for periodicity or 
quasi-periodicity. Overall, the variability patterns in all optical bands are very well correlated 
with no discernable time lag between bands (see 

In the context of our WEBT campaign, 3C 66A was also observed at near-infrared 
wavelengths in the J, H, and K bands with the AZT-24 1.1-m telescope at Campo Imperatore, 



" http : //nedwww. ipac . caltech. edu/ 



Table 2. Inter-instrumental correction offsets for optical data. Only those telescopes are 

listed for which non-zero offsets were required. 



Observatory B V R I 

ARIES +0.12 
Belogradchik - - - +0.03 

Rozhen +0.06 

Torino -0.12 
Valle d'Aosta - - - +0.035 
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Fig. 3. — Details of the R-band light curve (no extinction correction) during the outburst 
around August 16, 2007. Particularly remarkable is the rapid decline on JD 2454331 (August 
18) with a slope of ~ 0.016 mag/hr. 
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Fig. 4. — Details of the R-band light curve (no extinction correction) on Oct. 28 - 29, 2007 
(JD 2454402 - 2454403). The rise around JD 2454402.8 has a slope of ~ -0.008 mag/hr. 
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Italy. The primary data were analyzed using the same standard technique as the optical 
data (see above), including flat-field subtraction, extraction of instrumental magnitudes, 
calibration against comparison stars to obtain standard magnitudes, and de-reddening. The 
sampling was not dense enough to allow an improvement of the data quality by re-binning. 
Also, the IR observations covered only a period of ~ 1 month so that we can not draw any 
conclusions concerning correlated variability with optical or radio fiuxes. 

The resulting K-band light curve is included in Fig. [T], and a comparison of all three 
IR band light curves (J, H, K) is shown in Fig. O Within the limited time resolution and 
coverage, the JHK variability tracks the variability in the optical bands. For example, small 
BVRI fiares on JD 2454387 and JD 2454391 - 2454392 are closely matched by the JHK light 
curves. The JHK light curves exhibit variability on similar time scales as the optical bands, 
but with somewhat smaller amplitudes of A{JHK)<j0.3 mag. Rapid variability with rise 
and decay rates of A{JHK)>^0.1 mag/day can be seen, in particular around JD 2454387 - 
2454392. 



2.2. Radio observations 



At radio frequencies, the object was monitored using the Noto radio telescope at 43 GHz, 
the 22 m radio telescope (RT 22) of the Crimean Radio Observatory at 36 GHz, the Uni- 
versity of Michigan Radio Astronomy Observatory (UMRAO) 26 m telescope, at 4.8, 8, and 
14.5 GHz, the 14 m Metsahovi Radio Telescope of the Helsinki University of Technology, at 
37 GHz, and the Medicina radio telescope at 5, 8, and 22 GHz. 

The observations at 43 GHz have been performed using the 32-m radio telescope in 
Noto (Italy). The observations reported in this paper were carried out using the 43 GHz 
supereterodyne cooled receiver w ith a 400 MHz in stantaneous band. For details of the data 
analysis and fiux calibration, see lOtt et al.l (119941 ). 



The observations at 36 GHz were carried out at the 22 m Crimean Astrophysical Ob- 
servatory radio telescope (RT 22), using two similar Dicke switched radio meters. For details 
of the analysis and calibration of data from the RT 22, see I VolvachI (120061 ). At the UMRAO, 
the source was monitored in the course of the on-going long-term blazar r nonitoring p r ograrn . 
The data were analyzed following the standard procedure described in lAller et al.l (119851 ). 
As mentioned above, the sampling was rather poor, and some individual errors were rather 
large due to source confusion problems with 3C 66B. 

The 37 GHz observations were made with the 13.7 m diameter Metsahovi radio telescope, 
which is a radome enclosed paraboloid antenna situated in Finland (24° 23' 38" E, +60° 13' 
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Fig. 5. — Light curves in the near-infrared J, H, and K bands during the campaign, compared 
to the R-band hght curve. 
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Fig. 6. — Radio light curves of 3C 66A in 2007. 
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05"). The measurements were made with a 1 GHz-band dual beam receiver centered at 36.8 



GHz. A detailed description on the data reduction and analysis is given in iTerasranta et al. 



fll998h 



The radio data at 5, 8, and 22 GHz from the Medicina antenna were taken within a 
long-term monitoring program of ga mma-ray br i ght b lazars. For details on the analysis of 



data from the Medicina antenna see iBach et al.l (120071 ). 



The 36 and 37 GHz light curves from the Crimean Astrophysical Observatory and 
Metsahovi, respectively, are included in Fig. [1] and in the compilation of all radio light 
curves throughout 2007 in Fig. [61 They reveal small-amplitude (AF/F<^20 %) variability 
on time scales of ~ 10-20 days. The sampling is insufficient to allow any conclusions about 
possible correlations of the radio variability with optical or infrared variability. 

The most dramatic radio variability appears to occur at the highest observed radio fre- 
quency of 43 GHz at the Noto radio observatory (see Fig. [6]), possibly indicating erratic 
variability on a one-day time scale. However, we point out that the observed erratic vari- 
ability at higher radio frequencies (> 14.5 GHz) may, at least in part, be a consequence 



of int erstellar scintillation (for a more detailed discussion of this aspect, see iBottcher et al. 



20051 ). 



3. Optical spectral variability 



To investigate optical spectral variability observed during our campaign, we first com- 
puted snap-shot optical continuum spectra for several points in time throughout the cam- 
paign. A compilation of several such optical / near-IR sections of t he SEDs are plotted in 
Fig. [71 As in our earlier campaign on 3C 66A (IBottcher et al.ll2005l ). we find the (presum- 
ably) synchrotron peak of the SED in the optical regime, around (4 - 6) x lO^"' Hz. In order to 
investigate whether there is a systematic shift of the synchrotron peak with varying optical 
fiux level, we calculated B - R color indices for any pair of B and R magnitude measurements 
taken within 20 minutes of each other. Fig. [HI shows the resulting plot of the B - R color 
index vs. time, together with the R band light curve. This figure does not confirm the hint 
in our compilation of snap-shot optical/IR SEDs of a correlation between the brightness and 
the hardness of the source during our campaign. 

This is further illustrated in Fig. [9l showing the hardness-intensity diagram of B - R 
color index vs. R-band m agnitude. No trend o f B - R index with brightness can be found. 
In our previous campaign (IBottcher et al.ll2005l ). we had found a weak positive correlation of 
hardness with fiux for low fiux states, i?dered^l4.0 (in agreement with the positive correlation 
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Fig. 7. — Intrinsic (after extinction correction) optical - near-IR spectral energy distributions 
during 6 epochs throughout the campaign. 
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Fig. 8. — Intrinsic (after extinction correction) light curves of the R magnitudes and the 
intrinsic B - R color index of 3C 66A over the duration of the entire campaign. No obvious 
correlation can be identified. 
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Fig. 9. — Intrinsic optical hardness - intensity diagram for the complete data set over the 
entire duration of the campaign. No correlation of the B - R index with R-band magnitude 
can be identified. 
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identified by Vagnetti et al.|[2003l ). but no correlation at brighter flux state s. Similar results 



were also found for the low-frequency peaked BL Lac object BL Lacertae (jPapadakis et al. 



20071 ). Since during the entire 2007-2008 campaign the source remained in a bright state 



with i?dered^l4.0, our results are in perfect agreement with our previous flndings. 



In lBottcher et al.l (120051 ). we had also found hints of systematic spectral evolution during 
optical flares, with peaks of the B - R hardness occurring consistently a few days before 
peaks in the R-band light curve. Hardness - intensity diagrams of an individual flare showed 
a tentative hint of spectral hysteresis. During our 2007-2008 campaign, we did not flnd 
such indications. Fig. [10] shows the time-resolved hardness-intensity diagram and spectral 
evolution during the exceptionally bright outburst around Sept. 15, 2007 (~ JD 2454348 - 
2454370). The points are labeled by their respective JD, and the arrows indicate the time 
sequence. Our results do not conflrm the existence of spectral hysteresis in 3C 66A. 



4. Inter-band cross-correlations and time lags 

Our spectral variability analysis in the previous section already indicated the lack of 
systematic spectral variability. This would be consistent with a tight correlation at zero time 
lag between diff erent optical bands. We tested this hypothesis using the discrete correlation 



function fDCF. Edelson fc KrolikI 119881 ) 



Fig. [TT] shows a compilation of the DCFs of the B, V, and I bands with respect to R as 
a reference band. We calculated DCFs with various time lag bins and obtained consistent 
results, independent of our choice of the bin size. In order to test for possible time lags 
between the various light curves, we fltted the DCFs with an asymmetric Gaussian of the 
form 

Table 3. Best-flt parameters of an asymmetric Gaussian (see Eq. [1]) to the discrete 
correlation functions between various optical bands. 



Band (vs. R) 


B 


V 


I 


Fo 


0.703 


1.168 


0.776 


<Tr [d] 


2.01 ±0.51 


0.74 ± 0.08 


2.02 ±0.55 


(^d [d] 


1.88 ±0.20 


1.28 ±0.06 


2.30 ±0.27 


rp [d] 


0.13 ±0.19 


-0.01 ±0.05 


0.09 ±0.21 
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Fig. 10. — Intrinsic optical hardness - intensity diagram for the flare around Sept. 15, 2007. 
Individual points are labeled by JD - 2454000, and the arrows indicate the time sequence. 
No hints of spectral hysteresis or other systematic spectral evolution are visible. 
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Fig. 11. — Discrete correlation functions (DCF) of B, V, and I band vs. R band light curves. 
Different symbols indicate different time lag bin sizes, as indicated in the legend, indicating 
that our results are independent of our choice of a bin size. The curves indicate our best fit 
with an asymmetric Gaussian. 
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DCF{r) = Fo{-~Z. io.r<r, 



for T > Tp 



where Fq is the peak value of the DCF, Tp is the delay time scale at which the DCF peaks, 
and ar,d parametrize the Gaussian width of the DCF. The best-fit functions are plotted in 
Fig. [HI and the best-fit parameters are quoted in Tab. [31 These results confirm that there 
is no evidence for time lags among the optical bands in our campaign data. The rather low 
normalization of the resulting DCF fits of ~ 0.7 and ~ 0.8 in the B vs. R and I vs. R 
correlations might be a consequence of uneven sampling with several substantial data gaps 
in the B and I light curves (see Fig. [T]). 



5. Discussion 

The location of the synchrotron peak in the optical regime allows an estimate of a 
combination of the magneti c field B = Bg G and the Doppler factor D = 1{) Di. We follow 



the approach of Eq. (8) in iBottcher et al.l (120031 ). Here, the magnetic field energy density 



is parametrized as a fraction of the energy density in relativistic leptons. Furthermore, 
constraining the size of the emission region by the observed variability time scale, Tvar.obs = 
1 Th hr, through R = CTvar.obs D/ {1 + z), we find 



1/7 



= 8-4 ( :;;3 1 g (2) 



Ajf-iQ e| (1 + z) 

where 0^27 = 7.38 is the luminosity distance in units of 10^'' cm, /_io ~ 0.5 is the peak vFy 
flux in units of 10~^° ergs cm~^ s~^, and esy,-6 = Csy/lO"^ ~ 5 the dimensionless photon 
energy, e = /iz//(meC^), at the peak of the synchrotron component at i^sy ~ 5 x 10^^ Hz. p 
is the spectral index of the nonthermal electron population radiating above the synchrotron 
peak. With these values, Eq. ([2]) yields 

Be. = 19er(p-2)-VV-^/^Z^r^^/^G. (3) 

The lack of spectral variability may be interpreted as evidence that the synchrotron 
cooling time of electrons radiating in the optical regime is shorter than the light crossing 
time, which (including a factor D/{1 + z)) then determines the shortest observed variability 
time scale. If this were not the case, frequency-dependent radiative cooling effects would 
be expected to lead to observable time lags and systematic spectral variability. Hence, the 
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condition r'sy {1 + z)/ D < Tvar.obs may be translated into a Doppler-factor dependent lower 
limit on the magnetic field. Regarding electrons emitting in the R band, we find 



rsy,obs = 7.32 X 10^ 
which yields the magnetic-field limit 



l + z 
D 



1/2 



B, 



-3/2 



s < r. 



var.obs 



(4) 



B > 3.9D-^/^T~^^^ G. 

Combining this limit with the magnetic-field estimate of Eq. 
the Doppler factor: 



(5) 

we infer an upper limit on 



D <28 {p-2)-^'-^'t, 



3/32 ^-1/8 3/16 



(6) 



This restricts the Doppler factor to values typically found in the modeling of blazar-type 
quasars and low-frequency peaked BL Lac objects and excludes extreme values of F ~ D >,50, 
found in one-zone, leptonic synchrotron-self-Co mpton models for several high-frequency 
peaked BL Lac objects d etected at E>J. TeV (e.g., Begelman et al.l[2008 : Ghisellini fc Tavecchio 
20081 : iFinke et al.ll2008bl ). This aspect is particu larly exciting now, a fter the recent detection 
of 3C 66 A at > 100 GeV energies by VERITAS jSwordv et allboosi ). The very high Doppler 
factors mentioned above are primarily required to accomodate the very rapid variability, on 
time scales down to just a few minutes, at GeV - TeV energies seen in a few HBLs. Should 
future VHE 7-ray observations of 3C 66A indicate such rapid variability for this object as 
well, this would pose serious problems for conventional blazar emission models. 



6. Summciry 

We reported on an extensive optical, near-IR (JHK), and radio monitoring campaign 
by the WEBT throughout the fall and winter of 2007 - 2008, prompted by a high optical 
state in September 2007 with R < 13.5. Twenty-eight observatories in 12 countries in North 
America, Europe, and Asia contributed to this campaign. 

The source remained in a high optical state throughout the observing period and ex- 
hibited several bright fiares on time scales of ~ 10 days, including an exceptional outburst 
around September 15 - 20, 2007, reaching a peak brightness at i? ~ 13.4. Our campaign 
revealed microvariability with fiux changes up to |di?/dt| ~ 0.02 mag/hr with occasional 
indications of a slower rise than decay of individual fiares. 
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Our observations do not reveal evidence for systematic spectral variability in the overall 
high state (i?dered^l4.0) covered by our campaign, in agreement with previous results for 
the same overall brightness state. Note, however, that a positive correlation between optical 
spectral hardness (B - R color index) and optical flux at dereddened R-band magnitudes 
of i?dered^l4.0 has bccu fouud before. In particular, we do not find evidence for spectral 
hysteresis in 3C 66A for which hints were found in an earlier campaign in a somewhat lower 
flux state. On the same note, we did not find evidence for time lags between optical bands. 

During our campaign, as in earlier campaigns, the peak of the low-frequency (presum- 
ably synchrotron) component of the SED of 3C 66A was located in the optical regime. 
From the observed synchrotron peak flux, we infer a value of the magnetic field of the emis- 
sion region of -B ~ 19 e'^ D^^^'^ G. From the lack of systematic spectral variability, a 
Doppler-factor dependent lower limit on the magnetic field can be derived, which implies an 
upper limit on the Doppler factor, D < 28r^ ^^'^ . This excludes extreme values of the 
Doppler factor of i5>,50, inferred for some high-frequency-peaked BL Lac objects detected 
at TeV energies. 
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